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Abstract. We study the subsonic flows governed by full Euler 
equations in the half plane bounded below by a piecewise smooth 
curve asymptotically approaching xi-axis. Nonconstant conditions 
in the far field are prescribed to ensure the real Euler flows. The 
Euler system is reduced to a single elliptic equation for the stream 
function. The existence, uniqueness and asymptotic behaviors of 
the solutions for the reduced equation are established by Schauder 
fixed point argument and some delicate estimates. The existence of 
subsonic flows for the original Euler system is proved based on the 
results for the reduced equation, and their asymptotic behaviors 
in the far field arc also obtained. 



1. Introduction 

In this paper, we study subsonic polytropic flows governed by two- 
dimensional steady, full Euler equations: 

( V • m = 0, 

_ fxa (g) m\ _ 

V-f— — J+Vp = 0, (1.1) 

{ V-(m(£ + p/p)) = 0, 

where V is the gradient in x = (xi,^) G M 2 , m = (mi, 7722) the 
momentum, p the density, p the pressure, and 



E 



|m| 2 p 



2p2 ( 7 - i) p 

the energy with adiabatic exponent 7 > 1. The sonic speed of the flow 
is defined by 

c = V ipfp- 
The flow is said to be subsonic if |m/p| < c. 



2000 Mathematics Subject Classification. 35J15; 35M20; 76G25. 
Key words and phrases. Subsonic flows; full Euler equations; polytropic gas; 
nonlinear elliptic equations. 

1 



2 



JUN CHEN 



To my best knowledge, no theoretical result was obtained for sub- 
sonic flows governed by full Euler system in unbounded domain. There 
are rich literatures of subsonic potential flows, which is a simplified 
model for Euler flows. Shiffman obtained the first existence result in 
[31], using variational method. In [1], Bers used complex analysis to 
show existence and uniqueness for the subsonic potential flows. Finn 
and Gilbarg [T71 [18] solved the problem by PDE approach. Recently, 
Chen-Dafermos-Slemord-Wang [4] pushed the subsonic flows to the 
sonic limit, using the framework of compensated compactness. With 
in the same framework, Chen-Slemord-Wang [12j obtained transonic 
solutions by a vanishing viscosity method. Other results for subsonic 
or transonic flows of various models can be found in [2j [3], [5J [7J 

[3ra[IIl[l3l[H[23[20l[2aH provide 
related background and introduction. 

The domain we study is the upper half plane with piecewise smooth 
boundary asymptotically flattened as \xi\ — > oo. This setting can be 
viewed as a symmetric airfoil problem. For a airfoil symmetric about 
a;i-axis, we cut the exterior domain in half along the symmetry axis, 
and the upper part becomes our domain. More generally, we allow the 
boundary to be curved away from the profile. In this way, our setting 
also includes the model for the wind glancing the landscape. 

Let U = (m, p, p) be the solution for the subsonic flow. We prescribe 
an asymptotic limit Uoo, close to a constant subsonic state Uo, for 
the flow in the far field. Unlike the setting for potential flows, the 
asymptotic behavior is not a constant state. Otherwise, the full 
Euler system can be reduced to a potential flow (cf. Proposition 13. ip . 
To guarantee the convergence of the flows to in the far field, we 
need to obtain some decay property for if> — I, which is the difference 
between the stream function and its limit behavior. For the whole 
plane, one knows that the fundamental solution for a Laplace equation 
has the form log|x|. Therefore, in general, we can not expect the 
decay of a solution for an elliptic equation as |x| — ► oo in the exterior 
domain of the whole plane. This is the main technical obstacle for us to 
obtaining the subsonic flows in the whole plane. However, when flows 
are restricted in the half plane, we can exclude the logarithmic growth 
of solutions by prescribing proper decay condition at the infinity. 

We reduce the Euler system to a single elliptic equation (I3.19P for 
a stream function ip by capturing some conservation properties of the 
system. More precisely, three properties are contained in (13. 19ft : (1) 
existence of ip represents conservation of mass; (2) we use the fact 
p/ p 7 is constant along streamlines during the reduction, which implies 
entropy is conserved on each streamline. (3) to solve for p in terms 
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of if), V?/>, we use Bernoulli's law, which relates to the conservation of 
energy. Usually, stagnation points occur in various situations and cause 
major difficulties (for instance, regular reflection for Euler equations in 
self-similar coordinates). Our reduction process enables us to bypass 
the difficulty and to obtain the existence of solutions. However, we do 
not have uniqueness result due to the existence of stagnation points 
and complex behaviors of streamlines. More details are explained in 
Remark 16.11 

Once the Euler system is reduced to the elliptic equation (j3.19p . the 
remaining work is to solve this nonlinear equation. In detail, we first 
truncate the original domain Q by ball Br(0) centered at the origin 
with radius R. We solve (I3.19P in bounded domain Qr = Q D Br{0) 
with properly prescribed boundary condition. It is a standard method 
that we linearize (13.191) . construct a map T by solving the linearized 
equation, and prove the existence of a fixed point for T by Schauder 
fixed point theorem. The fixed point if>R is the solution for (13.191) in 
domain Qr. 

To take the limit of {iPr} as R — > oo and obtain the solution in Q, 
the estimates should be independent of the radius R. It makes the 
estimates complicated that there is no sign for the coefficient bo in the 
linear equation (14.11) . By choosing a proper barrier and using maximum 
principle with no restriction on the sign of bo (Lemma 14.21) . we obtain 
uniform estimates, independent of R, for ipR. The barrier function we 
construct only works for the half plane, not the whole plane. Whether 
one can find a suitable barrier function for the whole plane is unclear 
at this moment. 

The rest of the paper is organized as follows. In Section 2, we set up 
the subsonic flow problem, introduce the weighted norms, and state the 
main result. Section 3 explains the reduction of the full Euler system 
to a single elliptic equation for the stream function if). In Section 
4, we prove a technical lemma for the linearized equation and obtain 
crucial esitmates. In Section 5, we construct a iteration scheme to 
solve the nonlinear equation (13.191) in truncated domain Qr. Schauder 
fixed point argument is used to prove the existence of the solution. 
In Section 6, we take the limit of subsequence of solutions i/)r in Qr 
and obtain the solution in the half plane Q. The relation between the 
original Euler system and the reduced system is explained. 

2. Setup of the Subsonic Flow Problem 

In order to describe the conditions and results of our subsonic prob- 
lem, we need to use the following weighed Holder norms: for any 
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x, x' in a two-dimensional domain and for a subset P of <9fi, de- 
fine 5 X := min(dist(x, P), 1), <5 XX / := min(5 x , 5 X >, 1), A x := max(|x|, 1) 
and A xx / := max(|x|, |x'|, 1) . Let a G (0,1), a, /3 G R, and be a 
nonnegative integer. Let k = k 2 ) be a integer- valued vector, where 
hi, k 2 > 0, |k| = h + k 2 and L> k = d k x \d k x \. We define 

££kn = sup (^ (fe+CT ' 0) An^(x)l), 



L U Jfc,0;(/3);n 



x e o 

Ikl = fc 



r i(ct;P) / x max(fc+a+ CT ,0) A /3+fc+a \D k u(x.) - D k u(x.') 



x.x' e n 

x^x' 
|k| = k 



\\ U \\k,a;(p)-,(l — ^ [ M ]i,0;(/3) ; n + Mfc,a;(/3);n- (^-l) 
i=0 

For a vector-valued function u = (iti, w 2 , • • • , «n), we define 

n 

II U llfc,a;(/3);C / ., II II k.a:(3Y.Vf 
i=l 

Remark 2.1. In the definition of the weighted norms, the lower index 
in the parenthesis represents the weigh at the infinity and the upper 
index represents the weight to the set P, which will be the set of some 
corner points on the boundary in the paper. 

Define 

Ct-U Q ) = i u ■ IMlKkn < °°}- ( 2 - 2 ) 
For the weighted norms of functions in one-dimensional space T = 
(a, b), with either a = — oo or b = oo, we define 

[f] kMm = sup(|x| + l) fc ^|/( fc )(x)| 



> k+Q+ p\f {k) (x)-f ik) (x')\ 



[/]fe,a;(/3) ; r = sup (max( \x I , \x'\ ) + 1 

X,x'&,Xy^x' 

k 

||/IU,«iO0);r = 5^[/]fc,Q;(/3);r + [f]k,a;(J3);r ■ (2.3) 
i=0 

Our domain Q is the upper half plane bounded below by a piecewise 
smooth curve consisting of three parts: 

dn = T_UAUT + . (2.4) 

The following is the description of the three parts for dfl (see figure [1]) . 
Let r± = {x 2 = f±(xi)}, where /_ is defined on (—00, —1) and / + is 
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Streamlines, - - - " 




Figure 1. Domain for Subsonic Flows 

denned on (l,oo). Both T_ and T + approach xi-axis as |xi| tends to 
oo. More precisely, we let 

||/-||2,a;(a+/3);(-oo,-l) < 1) ( 2 -5) 
||/+l|2,a;(a+/3);(l,oo) < 1- (2-6) 

Let A_ = (-l,/_(-l)) and A + = (1,/+(1)) be the end points of 
T_,r + , respectively. The arch A connecting T± at A± can be param- 
eterized by 

f(s) = (f 1 (s),f 2 (s)), SG(-1,1), (2.7) 

where f(— 1) = A_, f (1) = A + and f\, fi are C 2,a smooth functions. 
We assume the angles 6*± between r ± and A at points A± satisfies: 

5<e° ± <7r-5, (2.8) 

for a fixed constant 5. 

Remark 2.2. The above condition guarantees that stream function ip 
for the flow is C 1,a up to the corner points, which means the flow U 
is C a up to the corners. If we allow corner angles 8± > n, ip will be 
C a up to the corners, and U will blow up at corners. We exclude the 
latter situation just to avoid unimportant details. 



Without loss of generality, we may also assume that 

^CB DO (0), (2.9) 
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i.e., A is contained in the ball of radius D centered at 0, 

f-M > ~\, f+M > ~\, /2(a) > -\- 

This means that domain fl is above the line xi = — |. 
We prescribe slip condition on the boundary: 

m.i/|an = 0, (2.11) 

where v is the outer normal on boundary dQ. 

Let (m*, 0,po, Po) be a constant subsonic solution for ( 11.11) . i.e., m*/p 
less than a/po/po- Fix constants po, Po and let mo < m* be a sufficiently 
small constant to be determined later. So Uq = (m , 0,po, Po), as our 
background state, is also a subsonic solution for (11. ip . We define a 
vector- valued function = (m^, 0,p , Poo) of variable x 2 as the as- 
ymptotic state for our solution U = (m, p, p) at the far field. We 
assume that is a small perturbation of the background solution Uq\ 

Woo — £^o||2,a;(0);(0,oo) < £™0i (2-12) 

where < e < \ is a small parameter to be determined later. 
Set P = {A-, A + } as the set for the weight. 

Now we state our main theorem about the existence of the subsonic 
flow in the half plane Q: 

Theorem 2.1. Suppose the boundary dQ satisfies (I2.5p -( I2TTU1) andUoo 
satisfies (12.121) . We fh ; < (3 < a < 1, depending on 5 in (12.81) . For 
sufficiently small m , depending on m*,p , p ,8,a, (3 and the profile 
A, and sufficiently small e, depending on m*,p , Po, ^ a , A A and m , 
there exists a subsonic solution U G C[ (Q) for (jl.ip with boundary 
condition (12.111) . such that 

\\U ~ U^;:^ < Cm , (2.13) 

where C is a constant only depending on m*,p , Po, 3, a, (3 and the pro- 
file A, but independent of m and e. 

Remark 2.3. Estimate (12.131) immediately gives the asymptotic behav- 
ior of the flow U. That is U approaches in C° norm at the rate 
|x| _/3 as |x| — > oo. 

3. Reduction of the Euler System 

In this section, we use the conservation properties of the Euler equa- 
tions (11.11) to reduce the four-equation system to one elliptic equation. 



and 

(2.10) 
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By the conservation of mass (first equation of ( 11. II) ). we can find a 
potential function ip for the vector field (-m2,mi), i. e., 

ijj Xl = -m 2 , ip X2 = m 1 . (3.1) 

From (11.11) . we can derive 

(mi9 xi +m 2 «9 ;C2 )(7lnp-lnp) = 0, (3.2) 

which implies that the quantity p 7 jp is constant along streamlines, 
provided that the solution is C 1 smooth. This constant only depends 
on the stream function ip . Thus, we have 

p = I^IaMp*. (3.3) 

7 

We will determine function A later by U^. 

From (11.11) . we can also derive the Bernoulli's law: 

|m| + , T ^ = B (3.4) 



2p2 ( 7 _ l) p 

along the streamlines, where B is the Bernoulli constant depending 
on ip. With equation (13.31) and (13.11) . the Bernoulli's law (13 .4p can be 
written as 

i|V^i 2 + ^(V)p 7+1 = S(^)p 2 , (3.5) 
In the subsonic region, we have 

|VVf < c 2 p 2 = (7 - l)A{^)p^\ (3.6) 
Inequality (13.61) and the Bernoulli's law ( 13.51) implies 

^ > TO (3 ' 7) 

Let x = ||V^| 2 and h(p,ip) = B(^)p 2 - A(^)p 7+1 - Therefore, in 
subsonic region, 

dh , ( 2B A 

— = pA[ - — - p 7 - 1 < 0. 

dp p V(7 + i)^ p ; 

Hence, we can uniquely solve 

X = h(p, = B^)p 2 - A(^)p 7+1 (3.8) 
for p = p(x, ip) by implicit function theorem. 



s 
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From (13.81) . we can easily calculate 



Px ~ ( 7 + l)ApT-25p' (3 ' 9) 
B' P - A' pi 

* = ( 7 + W-i- 2jB - ^ 



Therefore, we compute 

X\X 2 

(7 + 1)^P 7 - 2Bp 

Px 2 = Px( t Px 1 tpx 1 x 2 + i>x 2 1px 2 x 2 ) + P^x 2 

-^ Xl ^X lX2 - 1px 2 1Px 2 X2 + ^X 2 (B'p 2 - A'p^ 1 ) 



(3.11) 



(3.12) 



(7 + \)AfP - 2Bp 
Now we can reduce the Euler system into one equation. We replace 
m in the second equation of (II. ip with (—ip X2 , ip xi ) according to (13. ip . 
Multiplying the second equation of (11.11) by (^+l)Ap~ i — 2Bp, and using 
the expressions (13. lip and (13 . 1 2j) . we obtain the following equation: 

TpxMiM, WOVW., - WO) = 0, (3.13) 

where 

a u (V, W) = (7 - 1)^(^)P 7+1 - € 2 (3-14) 

a 12 {tp, V-0) = a 2 iO, Vip) = tp xi ipx 2 (3-15) 

a 22 (^, V^) = (7 - l)A^)p^ +1 - < (3.I6) 

FN, W) = ^^p 7+3 (7Afi' - 2A'B + AA'p^ 1 ). (3.17) 

7 

Similarly, the third equation of (11.11) gives rise to 

^(ciijN), VVOV^ - WO) = 0. (3.18) 

For a system without stationary points, i.e., Vip is nowhere 0, the 
original Euler system ( 11. ip can be reduced to the following equation 
for subsonic flows: 

ay(^VV0^ = ^(V>,V# (3.19) 
Equation (13.191) can be written in divergence form: 



V . (Yf) =B'p--A'p f . (3.20) 
\ P J 7 
Now we use the limit function to determine A, B and the limit 
function Z(x 2 ) of the stream function ip as |xi| — > 00. 
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Define 

1-X2 

Ifa) = / m 00 (s)ds. (3-21) 
Jo 

By (I27L2]) . we know that 

^m < = I' < 2m , 
which implies that I is invertible and 

^m x 2 < l(x 2 ) < 2m x 2 - (3.22) 



Let 



A( \ 7P0 



(7 - l)Poo(^2) 



b(x 2 ) = m B X2 l + 7Po 



Then define 



A(s) =A(r 1 (s)), S(s) = B{r\8)). (3.23) 



To describe the properties of A and -B, we need to modify the 
weighted norm in (I2.3P as follows: 

[/W);r = sup(|x|+mo)^|/«(x)| 

[f]' k , am r = snp (max(|xUx1) + mo) fc+ ^ l/(fc)( , X) J^ )(3:/)l 

ll/||fc,a;09);r = 5^[/]fc,O;09);r + Lf]&,a;(/3);r- (3.24) 
i=0 

Basically we replace 1 in the weight in (12.31) with m for the scaling 
reason. 
Set 

, IPo R _ m o 7Po 
A — 7 TT^y) -do — 7r~9 + 



(7-lK' 2pg ( 7 -l)Po' 

By (I2.12p . we conclude that 

\\ A ~ Aj|r 2 ,ai(o) 5 (o,oo) ^ C em , (3.25) 

II- 6 - £o|| 2 ,a;(0);(0,oo) < Co^ , (3.26) 

where Co is a constant depending only on m*,p ,p . 

Let us discuss the asymptotic behavior of U as |xi| — > oo. We do 
not expect constants states at the infinity for general subsonic flows 
governed by full Euler equations. Actually, if the flow is uniform at 
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the infinity, we only get a potential flow. This fact is described by the 
following proposition. 

Proposition 3.1. Suppose U is a C l solution of (11.11) with no stag- 
nation points fm is nowhere 0). If is a constant state, the flow U 
is potential, i.e., the velocity u = m/p is irrotational. 

Proof. If Uoo is constant, we immediately get A, B are constants by the 
procedure of obtaining A, B. Hence, we have A' = B' = . Equa- 
tion (13.201) . which is equivalent to (11.11) under the assumption in the 
proposition, becomes 




Since 

4> X1 = -m 2 = ~pu 2 , ip X2 =m 1 = pui, 
the above equation is just the irrotationality condition for the velocity 
( Ul )x 2 ~ ( U2 )xi = 0- Therefore, we have a potential flow. □ 

In general, the Euler system (II. ip and equation (I3.19P are not equiv- 
alent, because the streamlines may not be nice enough for us to do the 
reduction of the system. However, the solution of (13.191) guarantees 
the existence of the solution for Euler equations (11.11) . Therefore, we 
only need to solve (I3.19P in order to prove Theorem 12. 11 

By the definition of the stream function x/j, the slip condition (12.111) 
becomes the Dirichlet boundary condition for equation (13.191) : 

ip\aa = 0. (3.27) 
We define S set for the solutions of (I3.19P : 

Z = {u:\\u-l\\^-^ <C*m }, (3.28) 

where constant C*, depending on m*,p , p , a, (3, 5, A, will be deter- 
mined later in the estimates. 

We state the following theorem, which implies Theorem 12.11 

Theorem 3.1. For sufficiently small mo , depending onm*, p , p , a, (3, 6, A, 
and sufficiently small e, depending on m*,p , p , a, j3, 5, A andm , there 
exists a unique solution for equation (I3.19P with boundary condition 
(I3.27P in the set S defined in (13.281) . 

We split the proof of Theorem 13.11 into the following steps: 

(1) Use bounded domain Qpt := Q n Br{0) to approach Q, where 
Br{0) is the ball with radius R and centered at the origin. We 
linearize equation (I3.19P and solve the linear equation in Qr. 
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(2) With proper estimates for the linear equation, we solve the non- 
linear equation (13.191) in bounded domain fl R using Schauder 
fixed point theorem. 

(3) Let R — > oo, we prove the existence of solution for (13.19P in Q. 

(4) Estimate the difference of any two solutions in Qr and then let 
R — > oo to obtain the uniqueness of the solution for (I3.19p . 

Henceafter, we will use C to denote generic constants, depending on 
the fixed data m*,pQ, Po,a, /3,5, and the profile A, but independent of 
m , R. 

4. Estimates of a linear equation 

In this section, we study a linear elliptic equation 

aij(-x)u x . x , + bi(x)u Xi + b (x)u = 0, (4.1) 

in domain Qr. The estimates of this equation will be used later for the 
linearized equation. For equation (14.11) . we have the following assump- 
tions for the coefficients: 

",,(x)C<', > K^i + &) for any £ E R, (4.2) 

3 

\\aij - e6ij\\c<*(si R ) + < Cm , (4.3) 

i=l 

(x 2 + l)(|6i| + \b 2 \) + (x 2 + l) 2 \b \ < Cm , (4.4) 

where A, e are constants depending on p , p , and 6ij = 1 for i = j, 
otherwise Sij = 0. 

We let R > D + 1 so that the boundary of Qr includes the whole 
profile A. Let Sr = {|x| = R}C\dQR. Now we not only have the corner 
points A_,A + , but also have additional corners as the intersection of 
Sr with dQ. These two points are denoted by 

s_ — r_ n Sr, Si = r_(_ n Sr. 

Then we define the set of boundary points for the weight as: 

P = {A_, A + , S_, S+}. 
The boundary condition for (14.11) is: 

u\an R = g, (4.5) 

where 

IMI^-SU < Omo. (4.6) 
We then have the following lemma 
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Lemma 4.1. Suppose u G C 2,0: (VLr) fl C(VLr) is a solution for (14.11) 
with boundary condition (14.51) and assumptions (14. 2 p - (14. 6 p hold. For 

sufficiently small mo independent of R, we have the following estimate 
for u : 

IM&)£i < c *™°> ( 4 - 7 ) 

where C* is a constant independent of mo, R- 

To prove Lemma 14.11 we need a maximum principle for the elliptic 
equation (14. ip without restriciton on the sign for b^. We take the 
following lemma from [21] (Theorem 2.11): 

Lemma 4.2. Let elliptic operator L = aijd Xi d Xj + bid Xi + b$. For any 
bounded connected domain T>, assume a^^bi G C°(T>) and Oy satisfies 
ellipticity condition (14.21) . Suppose there exists a function v G C 2 (T>) fl 
C l (T>) such that v > in T> and Lv < in T>. Suppose u G C 2 (T>) fl 
C{V) satisfies Lu > inV. Then - achieves its nonnegative maximum 
on the boundary dT>. 

With this maximum principle, we start to prove Lemma [4.11 

Proof. Notice that there are two different weights in the norm in (14.71) 
(See definition of weighted norm (12.11) ) : the weight with upper index 
—a — 1 is for the small scale near the corner points A_, A + , S^, S+ 
and the weight with lower index (3 is for the large scale away from 
the origin. We split the proof into three parts: Part 1 is for the 
estimate of maximum norm of u in the whole domain Qr, Part 2 is 
for the region near A; and Part 3 is for the region far away from the 
profile A. Let D = 2Dq + 1, where Dq is the radius to bound the 
profile A. Let Qd = H {|x| < D} be the region for Part 2, and 
Vt c D = Q R n {|x| > D - 1} for Part 3. 

Part 1. In this part, we first construct a comparison function v and 
use the maximum principle (Lemma 14.21) in the whole domain to 
obtain the control of the maximum norm of u. 

Define the comparison function v by 

v (x) = r~ a - (x 2 + l) a , (4.8) 

where r = \Jx\ + (x 2 + l) 2 . 

We now verify the fact that Lv < 0. 
First, it is easy to compute 

Av = ((3 2 - a 2 )r- a - f3 - 2 (x 2 + l) a 
-a(l - a)r- a -P{x 2 + l) a ~ 2 
< -a(l - a)r- a - p {x 2 + If" 2 , (4.9) 
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noticing that < /3 < a < 1. Also, one can verify that 

\Dv\ < Cr- a -\x 2 + l)"" 1 , (4.10) 
\D 2 v\ < Cr- a -\x 2 + l) a ~ 2 . (4.11) 

We rewrite Lv as 

Lv = (L- eA)v + eAv. 

By assumptions fl4.3[) and (14.41) . together with (I4.10p and (14.111) . we 

have 

\(L - eA)v\ < Cm r~ a -P(x 2 + 1) Q " 2 . 

The above estimate and (14. 9 p imply that Lv < in Q R , provided m 
is small enough. Obviously, v is positive. Hence, by the maximum 
principle (Lemma 14.21) . we conclude that 

u ^ \g\ 

— < max — . 

v dn R v 

By replacing v with — v and using Lemma 14.21 again, we have 

\u/v\ < max \g/v\. 

dtt R 

This, with assumption (14.61) . implies 

|w(x)| < Cm r~ p . (4.12) 

Part 2. For the region near the profile A, we need to take care of 
the corner points A_,A + . We use the weight up to P and drop the 
lower index j3 for the weight away from A. We treat the corner v4_ 
first, and A + can be dealt in the same way. For convenience, we move 
A- to the origin O. Assume the angle between T_ and Xi-axis at O 
(original AS) is 8-, and the angle between A and Xi-axis at O is 8q. 
Therefore, the tangential directions of T_ and A at O are 

z/_ = (cos 6L, sin 6L), v = (cos 6*0, sin 9 ), 

respectively. Let u = u — g{0) — C\x x — c 2 x 2 , where Ci,c 2 are linear 
combinations of ^-(O) and §^{0) through solving the linear system 

r ( Ci , C2 )-^ = !ho) 

\ ( Cl ,c 2 )-u = ^(O). 

Hence 

|ci| + M <C\Dg(0)\, 

and 



u(O) = 0, Du(0) = (0,0). 



(4.13) 
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Choose r > small enough, such that f2£>n£> ro (0) is connected. Thus, 
for this fixed radius r < min(_D ; 1)> 

u(x)\ d n DnBro (o) < Cm |x| 1+a . 

We know u satisfies the following equation 

Lu = ciijdidjU + bidiU = F , (4-14) 

where F = —biCi — b u. By estimate (14. 12ft and condition (14.41) . we 
have 

|F | < Cm . 

Notice the elliptic operator in (I4.14p does not contain bolt term. So 
we can use standard maximum principle to control u. The comparison 
function for u is defined in polar coordinates (r, 9) by 

Vl (r, 6) = Cm r 1+a sin(r + (6 - 9 )), 
for small positive r depending on 6- — 8q and a. One can check that 

Lvi < —CrriQ < F = Lu. 
Also the boundary condition satisfies 

vi\d(n D r\B r0 (O)) > Cm r 1+a > u\d(n D nB ro (o))- 
By maximum principle, we conclude that 

\u\ < Cm r 1+a } 

for |x| < r . 

Once we have the above estimate near the corner A_, we use Schauder 
estimates with proper scaling to obtain the estimate near the corner 
A-{0): 

ll^lU^U) < Cmo. (4.15) 

The procedure is standard and related details can be found in chapter 
6 of |22j. One can also refer to [5] (Lemma 4.2) for similar scaling 
argument. We sketch the proof as follows. 

For any x G Qd H B ro / 2 , let the angle between T_ and the ray A_x 
be 9 X0 , and the angle between T_ and A be 9°_. Consider two cases: 
Case 1, # XQ > 7r/6 and 9- — # Xo > 7r/6 ; Case 2, otherwise. For Case 
1, we know that u satisfies equation 

Lu = -hd - b (g(O) + cixi + c 2 x 2 ). 

Take the ball S| Xo |/ 2 (x ) C ftp fl as the domain and by Schauder 
interior estimate (Theorem 6.2, [22]), we have 

ll^ll2^ |xol (xo)<^o|x o r +1 . (4.16) 

— T~ 
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Here the upper index (0) is understood as the weight up to dB\* \ (x ). 

2 

For Case 2, let x* be a boundary point with the shortest distance 
from x , and d* = |x | sin(|#°J. Hence, i?7d^(x*) still contains x*. 

8 

We use Schauder boundary estimate (Lemma 6.4, [22]) in the domain 
Qd n Bd*(x*) to obtain the estimate: 

NI^M < Cmolxor 1 . (4.17) 

Combining (14.161) and (14.171) gives the corner estimate (14. 15j) . The 
other corner A + is treated in the same way. Together with standard 
Schauder estimates away from the corners, we conclude the estimate 
in Qjj: 

Nl£S P) < Cm o- (4.18) 

Part 3. For the domain Q C D , we also consider two kinds of estimates: 
one is near the corner points S 1 ^, S+, the other is away from the corners. 

The corner estimates are similar to those in Part 2. In brief, consider 
the corner for instance. If x £ Q c D r\B rQ (S^), we have the following 
the estimate 

ni£;j^W) ^ Cm oR- . (4.i9) 

If xo G £l c D and |xo| < R/2 , the ball f?|x i (xp) has no intersection 

2 

with the outer boundary Sr = {|x| = R} or the profile A. Using 
conditions (I4.2IUI4.3I) . (14.61) and estimate (I4.12p . by Schauder interior 
estimates (see Theorem 6.2 in 1221), we have 



ll«ll2,o;08);B|^|(xo) < Cm ^ ( 4 - 20 ) 

where no upper index in the norm means no weight up to P. For 
l x o| > -R/2, we use Schauder boundary estimates (Lemma 6.4 in [22] ) 
with the boundary condition H4.6[) and estimate (14.191) to obtain 

Estimates (14.181) . (14.201) and (14.211) imply estimate (1477]) in the lemma. 

□ 

By the continuity method, one can prove the existence of solutions 
for (14. ip with estimate (14. 7p . The uniqueness is simply the result of 
the maximum principle, Lemma 14.21 with the aid of the comparison 
function v constructed in Lemma |4~T1 Since the procedure is standard, 
we omit the proof and only state the result as follows. 
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Lemma 4.3. Assume (14.2p — (14.60 holds. For sufficiently small mo, 
equation (14. lj) with boundary condition (I4.5|) admits a unique solution 

5. Nonlinear Equation in Bounded Domain 

In this section, we will solve equation (13.191) in the bounded domain 
Qr with boundary condition given below. 

We want to prescribe the boundary data for ip such that (13.271) hold 
on dQ fl and ip — I vanishes away from dfl. We will define function 
g such that 

ip — I = g on d£l R . (5.1) 

First it is easy to construct a smooth cutoff function r](s) such that 
rj(s) = 1 for |s| < Dq and r](s) = for \s\ > D + 1. We also assume 
that 1 1 r] | |c 2 .«(R) < 10. Let 

#(x) = -r)(x 2 ) ((1 - vM) U/sign^o^i)) + Z ( x 2)) , (5.2) 

where sign(xi) = — for X\ < and signal) = + for X\ > 0. 

It is easy to check that g\gn = — 1\qq, g = for x 2 > D + 1, and also 
g satisfies condition (14.61) . Now, let constant C* in (13.281) be the same 
as in estimate (14. T[) in Lemma 14.11 The set for solutions of (13. 19[) in 
Vt R is 

S R = {u : ||u - 1\\ { 2^)S R ^ C* m o}- (5.3) 
We state our lemma for the solution of (13. 19j) in VLr. 

Lemma 5.1. Equation (13. 19[) in VLr with boundary condition (15.11) 
admits a unique solution in the set provided m,Q,e are sufficiently 
small. 

Proof. To prove the lemma, we first linearize the nonlinear equation 
(I3.19p . By solving the linearized equation, we construct a map T in 
the set . The solution of the nonlinear equation (13.191) is a fixed 
point of T. 

We start the proof with the linearization of (I3.19p . We know that 
the limit function I satisfies (13.191) by its definition. That means the 
following equation holds: 

<kj(W)lxv*i = F(l,0,l'), (5-4) 
where ay, F are defined in (I3.14I) - (I3.16I) . Taking the difference of equa- 
tions (137191) and ((574]) leads to 

a -(V>, VV0(V> - 0*^ + l"Mip, W) - a 22 (i, 0, 1')) 
= F(^,V^)-F(/,0,r). (5.5) 
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Denote (/ + s(ip - I), V(7 + s(V> - 0)) b Y tjf , and let 

= Oy ( 



a« = a ii (V' ) V'0) ) (5.6) 



= i" l\a 22 )^ti)ds- I F^ds, 
Jo Jo 



(5.7) 



(6f,6j) = I" (a 22 )^(tf)ds- F^(tf)ds. (5.8 



Then we linearize equation (15.51) as follows: 

4$ - l) XiXj +bf$-t) + bi(4> -0 = 0. (5.9) 

We solve the above equation by applying Lemma I4~T1 In the following, 
we will check the conditions for ip e E R . In fact, ( Oil will 

be satisfied if (14.31) holds with sufficiently small mo- We let e in (14.31) 
be au(m x 2 , 0, m ) = 022(^0^2, 0, m ) = 7PoPo- By the expressions for 
<kj, fl334l) - fl3TT6l ). and (jOljl . fl3T9l) . it is not hard to verify that 

\\( aij Utt)\\c«(a R ) < — ccr, 

W(aij)x7ip(tf)\\c*(n R ) < — CC*, 

mo 

for ip E £ and mo small. 

In the proof of this lemma, the generic constants C are independent 
of C*. 

Let £ < , we have 

< / IIK-)^(tf)||c«(n fl )^||^ -m x 2 || c -(n fl ) 
■/ 

1 

1 1 (a<j ) (t^) 1 1 c-(n«)ds 1 1 V(^ - m x 2 )\\c<*(n R ) 

< eC{C*) 2 

< Cm . 

In the same manner, we can obtain 

II &i || cr* (0*) ^ Cm o, i = °> M- 
The above estimates lead to condition (14.31) . Now we verify condition 
(14.41) . For 6q , by its expression, we need to estimate {a 22 )^ and F^. By 
the definition of F, (13.171) . and estimates (13.251) . (I3.26p . we have 

em C 



\F^,Vij)\ < 



m + 
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Let u s = I + s(if) — I). For any if) G £, we consider two cases: Case 1, 
\ip — 1\ < ^m x 2 ; Case 2, otherwise. For Case 1, 

, 1 1 

u s >l- -m x 2 > -m x 2 , 

noticing ( 13.221) . Therefore, 

i ^ / ,/.m i ^ / ^ m em C eCm 

for e < mo 2 . For Case 2, i.e., \tf> — 1\ > ^m x 2 , since if> G £r, we have 

-m x 2 < \ip — l\ < 



x 



P 



This implies that x 2 < (AC*) 1 +p = R . Hence, for e < ((m /(l + R )) 
we have 

f ^ eC ^ eCm {l + Rq) 2 ^ Cm 



mo (m ) 2 (l + x 2 ) 2 (1 + x 2 ) 2 
The above analysis about both Casel and Case 2 gives rise to 

ml Fi(ti)ds - Cm ° 



< 7T^r (5-10) 



1+X2 



Similarly, we have 



I I < 



•i'2 



Together with the fact that \l"(x 2 )\ < j^r, we conclude that 



On < 



01 (l+x 2 ) 2 ' 

Same arguments apply to the estimates for bf,b 2 : 

Cm,Q 



11+ °£ <T— ■ 
1 + 2 2 

Therefore, we have verified condition ( 14.4ft . Thus, we apply Lemma 
14.11 to solve (15. 9p for ?/>, where u = if) — I in Lemma [4.11 By estimate 
( 14.71) . we know that ^ £ Sj. Therefore, we can define a map T from 
to itself by = if). It is obvious that a solution for the nonlinear 
equation (13 . 1 9[) is a fixed point of T. In order to prove the existence of 
a fixed point of T, we apply Schauder fixed point theorem, which says: 
if Y,r is a compact convex set of a Banach space £>, and T : — > 
is continuous in B, T has a fixed point. 
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Now we let B = C { ~^ P \n R ) (cf. definition (Q), where < a' < 
a. Obviously, E# is compact and convex in B. We only need to verify 
T is continuous in E^j. We prove this by contradiction argument. 

Suppose T is not continuous. Then there exists a sequence {if> n } C 
E R such that if) n — > if) in B, but Tip n does not converge to if) = Tip. 
This implies that we can find a subsequence {Tip nk } such that 

\\Tif) nk -if>\\ B >ca>0, (5.11) 

where || • \\g denotes the weighed norm || ■ ll^/^jyni anc ^ c o * s a nxe d 
constant. Since {Tifi nk } is compact in B, there exists a subsequence, 
still denoted by {Tifi nk }, convergent to if) e E#. On the other hand, 

ip n ^ ip in B implies that a i " k — > afj(i,j = 1, 2), 6- " fc — > &f (i = 0, 1, 2) 
in C a ' norm. Let k — > oo and we see that ^ is also a solution of (15.91) 
with the same boundary condition (I4.5p . Inequality (15. lip implies that 
— i>\\B > Co > 0, which means if), if) are two distinct solutions for 
(15.91) . This contradicts with the uniqueness of the solution for (15. 9|) . 
Hence, we verified the continuity of T in E#. 

By Schauder fixed point theorem, there exists a fixed point ipn for 
T. So ipa is a solution for the nonlinear equation (I3.19I) in Qr with 
boundary condition (I5.ip . 

The uniqueness of the solution for (I3.19P is proved by the maximum 
principle Lemma I4.2I as follows. 

For any two solutions if>i,if> 2 £ E fi of (I3.19p . (I5.ip . we take the dif- 
ference of the two equations and obtain 

a$kk- ^2% iXj + (^2)^(4-41) = Wu Wi) -F{ip 2 , Vif) 2 ). 
Similar to the notations in (I5.5l) - (l5.8p . we set 

t, = (if) 2 + s(if)! - if) 2 ), V(^ 2 + s(^i - ^2))), 
and let u = if)\ — if) 2 . Then we derive the following equation: 

+ hu Xl + b u = 0, (5.12) 



where 



bo = (^x^ / (a ij ) i ,(t s )ds - / F^(t s )ds, (5.13) 
Jo Jo 

(bi,b 2 ) = (fa)^ / (ciij) v ^(t s )ds - / F Vil (t s )ds. (5.14) 
Jo Jo 

Notice that the factor {i}) 2 ) XiXj in (I5.13P and (I5.14p blows up at the 

corner points in P = {A_, A + , S 1 ^, S 1 ^}. Lemma l4~2l requires continuity 

of coefficients a^,^ up to the boundary. Therefore, in order to apply 
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the lemma, we truncate small neighborhoods around corners from Q R . 
Define 

Bv,r = |J B r (I), (5.15) 
lev 

where V is a set of boundary points on Q R . Let 

ti% r = n R \B^:, s P!r = dB v>r nn R . (5.16) 



Now, we know 6j G C ^fi^' r J , i = 0,1,2. Define 77 = |x — J|, for any 
corner point I € P. Then we have the following estimates for bf. 

(x 2 + 1)(N + N) + (x 2 + l)>o| < C7m (l + J2 r i~^- ( 5 - 17 ) 

We construct a comparison function 5 as follows: Let 

Vl = rj~ 4 \x 2 + l)z. (5.18) 
Define v = YliieP v i- By (I4.9l) - fl4.11l) . we can verify that 

atHvj)^ < -c (rf^-\x 2 + 1)1 + rf^(x 2 + l)^" 2 ), (5.19) 



i=l,2 



< Cm (l + r( i' l ) r i |/? ( x 2 + l)^" 2 , (5.20) 

where c > is a constant only depending on (3. Inequalities (15.191) 
and (I5.20p directly imply that 

af!v x%Xj + biV x% + b v < 0. 



By Lemma [4.21 we have 



\u\ \u\ . . 

sup — = max — . (5.21J 



hai 



Since u\qq, r = and u G C ,a (iln), we know that 



max 

an£ r v s P,r v 



Hence, we have 



sup \u\ < CR~^r a sup v < CBJ p r a {l + r~s 



,P,r n P,r 



R 



Therefore, by letting r — ► 0, the above inequality implies sup Q |w| < 0. 
This shows the uniqueness of the solution for (13.191) in Q R . Hence, the 
proof of this lemma is complete. 

□ 
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6. Subsonic Flow in Half Plane 

After we solve (I3.19P in Qr, we let R tend to infinity to prove The- 
orem 13.11 as follows. 

Proof of Theorem IJ.il By Lemma 15.11 for a given radius R, we can 
find a unique solution ipR 6 E#. By a diagonal process, one can choose 
proper sequence R n — > oo as n — > oo, such that converges to some 
function -0 in || * H2 c/-(/3)'nQ norm > f° r an y fixed Q > D + 1. Since 
G we have 

ii^ n -C;S£L ^* m °' ((U) 

for any i? n > 2(D + 1). Let n — > 00 in (16.11) . we obtain estimate 

11^ " < C*m , (6.2) 

which implies ^ G S. This completes the existence of solutions for 

(HEED. 

To prove the uniqueness of the solution, we will use the asymptotic 
behavior of solutions described by set S. We still use the truncated 
domain Qr and follow the same strategy as in the uniqueness part of 
Lemma 15.11 Now, the situation here is slightly different from that in 
Lemma [5.11 (1) we have no singularity for (ip2)x i x :j at corners S^,S+; 
(2) u = ip\ — ip2 does not vanish on the boundary portion Sr = {|x| = 

R} nfi. 

Similarly as in Lemma I5.1[ we have the same equation (15.121) for u. 
We define the comparison function v by v = va_ + va + , where v a_ , va + 
are defined in (15.181) . The domain we consider here is fl^ r defined by 
(I5.15p . By the same computation as in Lemma [5.11 we have 

atjVxiXj + Ms* + b v < 0. 
We know that \u\ < CR'^ on Sr by the definition of S. Also, we see 

\u\s P J < Cr a , v\ Sr > Hrfr, v\ Sp<r > r-&. 
Therefore, we conclude that 

\u\ ^ 1 

sup^ < Cm&x(r a+ -^,R—^). (6.3) 

Let r = R~^ and we obtain 

sup|w| < CR'-^ sup v < CR-^(1 + Ri) < CR~i. 

f2 f2 
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Letting R — > oo gives rise to sup^ \u\ < 0, which implies the uniqueness 
of the solution for (13.1 9p in S. This finishes the proof of Theorem 
O □ 

Once we proved Theorem 13.11 define U = (m,p,p) by 

7-1 

m= (ipx 2 ,-'4 ) x 1 ), p = p(ip,Vip), p= A{ip)p~<, 

7 

where p is uniquely solved from Bernoulli's law (13.51) . Hence, by H3. 13j) . 
(I3.18P and (13.51) . we can recover the original Euler equations (11.11) . It 
is easy to check that U satisfies (12.131) with the aid of estimate (16. 2\i . 

Remark 6.1. The uniqueness of the solution for Euler system (II. ip can 
not be obtained from Theorem 13.11 due to two obstacles. One is the 
existence of stagnation points, which disqualifies equivalence between 
equation (I3.19P and the two momentum equations f)3. 13j) . (13.181) . The 
corners A_,A + on dfl are necessarily stagnation points, because Vif) 
is continuous up to the corners. Whether or what kind of stagnation 
points may appear inside domain Q is not clear. The other problem is 
about the complexity of streamline topology. During the reduction in 
section [3j we assume streamlines have simple topology, which means 
that streamlines in Q extend from — oo to oo in Xi, so that information 
about A, B can be carried along streamlines and reach the whole do- 
main Q. However, the geometry of profile A may be complicated and 
cause nontrivial topology of streamlines, such as closed orbits or inter- 
section of streamlines at stagnation points. The above reasons prevent 
us from obtaining the uniqueness for the Euler flows out of Theorem 
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